Simulation results for an interacting pair of resistively shunted Josephson junctions 
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Using a new cluster Monte Carlo algorithm, we study the phase diagram and critical properties 
of an interacting pair of resistively shunted Josephson junctions. This system models tunneling 
between two electrodes through a small superconducting grain, and is described by a double sine- 
Gordon model^. In accordance with theoretical predictions, we observe three different phases and 
crossover effects arising from an intermediate coupling fixed point. On the superconductor-to-metal 
phase boundary, the observed critical behavior is within error-bars the same as in a single junction, 
with identical values of the critical resistance and a correlation function exponent which depends 
only on the strength of the Josephson coupling. We explain these critical properties on the basis of 
a renormalization group (RG) calculation. In addition, we propose an alternative new mean-field 
theory for this transition, which correctly predicts the location of the phase boundary at intermediate 
Josephson coupling strength. 



I. INTRODUCTION 

The effects of dissipation and decoherence are ubiqui- 
tous in quantum systems and influence the properties 
of materials and nano-scale devices in a profound way. 
Already the simplest model system, an Ising spin in a 
transverse field which is coupled to an Ohmic heat bath, 
displays interesting behavior such as a dynamical phase 
transition to a localized state at a critical value of the 
dissipation strength^. Another prominent example is the 
resistively shunted Josephson junction, which undergoes 
a superconductor-to-metal transition at a critical value 
of the shunt resistance'^^ , which equals the quantum 
of resistance Rq = h/ie^ = G.Skfi. Arrays of Joseph- 
son junctions with dissipation have been studied both 
as a model for granular superconducting films or nano- 
wires^^^, and in their own right^"'^^. The behavior of 
these systems, and in particular their superconductor- 
to-metal phase transition, is far from being completely 
understood. 

Recently, Refael et al. have studied a model of a meso- 
scopic superconducting grain connected to two leads via 
Josephson tunneling and shunt resistors^. The phase 
diagram for this simple system was shown to be re- 
markably complex, with three distinct phases. In ad- 
dition, contrary to the case of a single resistively shunted 
Josephson junction, one of the phase boundaries is con- 
trolled in part by an intermediate coupling fixed point, 
and the superconductor-to-normal transition across this 
boundary can be tuned by the Josephson energy itself. 
The above effects indicate that although the system dis- 
cussed in Ref. 1 is zero dimensional, it nearly has the full 
complexity of a one-dimensional array of Josephson junc- 
tions. Much of the recent results on Josephson junction 
chains in Ref. 12 draw directly from the two-junction sys- 
tem. The RG-flow equations of the two-junction system 
are also nearly identical to those of the two-dimensional 
triangular lattice presented in Ref. 13. In this closely 
related work on Josephson junction arrays, the local na- 
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FIG. 1: Upper figure; Two-junction model considered in 
Ref. 1 with identical shunt resistors R for the left and right 
junction. The dotted box represents the central grain which 
incorporates a phenomenological charge relaxation mecha- 
nism described by the resistance r. 

Lower figure: Equivalent model with modified shunt resis- 
tors -Rjunction = R + 2r aud an additional resistor -Rioad = 
{R/r){R + 2r) connecting the leads. 



turc of the phase transition, as well as floating phases 
have been discussed. 

The development of a powerful new cluster Monte 
Carlo algorithm^** has allowed to test and verify seve- 
ral analytical predictions for the single junction and to 
observe continuously varying correlation exponents along 
the phase boundary. In this paper we will use adapted 
versions of these cluster moves to simulate the two- 
junction model of Ref. 1, which is the simplest exten- 
sion of the single junction case which exhibits interest- 
ing new physics. The model (with identical shunt re- 
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FIG. 2: Phase diagram in the hmits Ej / Ec -C 1 (solid hnes) 
and Ej/Ec ^ 1 (sohds and dashed hnes). Besides the su- 
perconducting (FSC) and metallic (NOR) phase, the two- 
junction system can be in a state which is superconducting 
from lead to load, although the individual junctions arc insu- 
lating (SC*). The parameters u, w and u, w are defined in 
Eqs. (22) and (26). The values m = and w = define the 
FSC-NOR boundary in the weak coupling limit, and w = 
and W = in the strong coupling limit. 

sistors) is shown in the upper part of Fig. 1. It con- 
sists of two Josephson junctions with coupUng energy 
Ej, each shunted with an Ohmic resistance R. On the 
central grain, the model incorporates a "charge relax- 
ation mechanism" (supposed to represent the break-up 
of Cooper pairs into electrons) which is described by an 
additional resistance r. 

Dissipation produced by electrons flowing through the 
resistors reduces phase fluctuations between the super- 
conducting islands which they connect. This can be seen 
from the dissipative action term in Eq. (4), and it ex- 
plains how strong dissipation leads to superconducting 
phase coherence. Depending on the values of R and r, 
three different phases occur: 

• The individual junctions are insulating and there is 
no super-current from lead to lead: Normal phase 
(NOR), 

• The junctions arc superconducting and thus also 
the whole device from lead to lead: Fully Super- 
conducting phase (FSC), 

• The individual junctions are insulating but there is 
superconducting phase coherence from lead to lead: 
SC* phase. 

The phase diagram for the limiting cases Ej <C Ec 
and Ej » Ec has been computed in Ref. 1 using an 

RG approach and is shown in Fig. 2 as a function of the 
resistances R and r. The boundary between the NOR- 



and FSC-phase in the region marked ICFP depends on 
the value of Ej/ Ec as well as the values of the resistors r 
and R. The behavior of the system in this region is con- 
trolled by an intermediate coupling fixed point (ICFP), 
which will be discussed in Sec. II D. 

In the remains of this paper we present the results of 
a thorough Monte Carlo (MC) investigation of the two- 
junction system. Each result is compared with predic- 
tions or explanations based on the RG flow equations pre- 
sented in Sec. II. After discussing the various phases of 
the two-junction system, we explain in Sec. Ill the Monte 
Carlo method which was used to investigate the model 
numerically. In particular, we discuss several types of 
eflicient cluster updates, which are adaptations of the re- 
cently developed rejection-free cluster algorithm for sin- 
gle resistively shunted j unctions. 

In Sec. IV we identify the three phases NOR, FSC and 
SC* by computing the temperature dependence of the 
lead-to-lcad and lead-to-grain resistance. The numeri- 
cally obtained phase diagram at intermediate Josephson 
coupling is compared to the theoretical predictions for 
weak and strong Josephson coupling. There is a good 
agreement between theory and simulation results, except 
in the region, where the three phases meet. We explain 
these small deviations of the measured phase boundaries 
from the predicted ones in terms of slow cross-overs in 
the RG flow, which prevent the MC simulation from ex- 
ploring the zero-temperature behavior. 

In Sec. V we concentrate on the critical FSC-NOR line 
in the ICFP region. In particular, we compare the critical 
properties of the two-junction system to those of a single 
jimction with the same Josephson coupling Ej/Ec, and 
a shunt resistance Rs — Rq. The two systems exhibit 
(within error-bars) identical behavior in their correlation 
functions, fluctuations and the effective resistance of the 
critical junctions. These critical properties depend on the 
value of Ej / Ec but not on the value of the resistors r and 
R. In Sec. V B, we try to explain this rather surprising 
observation in terms of the effective junction resistances 
of the two-junction system, calculated from the RG flow 
of Sec. II. We show that the predicted critical resistance 
exhibits only a weak dependence on the (shunt) resistors 
and agrees quite well with the measiired value. 

In Sec. VI, however, we pursue a different path to ex- 
plain the remarkable resemblance of the single resistively 
shunted junction and the two-junction system. We show 
that the data can be well explained by a 'mean-field' 
theory, which treats one of the two junctions as an ef- 
fective resistor. This is a new way of approaching the 
double sine-Gordon model action (Eq. 1 below). The as- 
sumption is that on the FSC-NOR phase boundary, each 
junction sees an environment which imitates a shunt re- 
sistor Rg = Rq and that it may be replaced by a resis- 
tor whose value can be found from the critical resistance 
of the single resistively shunted junction (with the same 
Josephson coupling). On the basis of these assumptions 
it is possible to derive an expression for the position of 
the FSC-NOR phase boundary, which fits the MC data 
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and the RG-based predictions quite well. 

In Sec. VI we consider the phase-phase correlation func- 
tions in the NOR phase within the ICFP region. For fixed 
resistors, we measure a strong dependence of the cor- 
relation exponents on the Josephson coupling strength, 
which is in contrast to the single junction model, where 
these exponents depend only on the value of the shunt re- 
sistor. This behavior is explained as a consequence of the 
flow in the additional Josephson coupling J-|_ between the 
leads, which is generated under the RG and not present 
in the single junction model. 

II. THEORY FOR THE SYMMETRIC 
TWO-JUNCTION SYSTEM 
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FIG. 3: In the weak-coupling limit, we consider the RG flow 
of the Josephson couplings J (with Jo = Ej) and J+. J 
is the amplitude for Cooper-pair hopping between cither of 
the leads and the grain. J+ is the amplitude for Cooper-pair 
tunneling from lead to lead. This process is generated in the 
second order of the RG flow. 



In this section we will present the cfFectivc action for 
the symmetric two-junction system, and then briefly dis- 
cuss its RG flow equations and the various phases. This 
discussion will prove to be especially useful when inter- 
preting the Monte Carlo results. 

We will first describe the NOR-SC*and SC*-FSC tran- 
sitions in the weak and strong coupling limits. It is im- 
portant to note that the SC*phase appears due to inter- 
actions between the two junctions, and can not be under- 
stood in terms of the physics of a single junction. Finally, 
we will discuss the important intermediate coupling fixed 
point which controls the direct NOR-FSC transition. 

A. Effective action 



The system discussed in Ref. 1 is equivalent to the one 
illustrated in the lower part of Fig. 1, where each junction 
is shunted by a resistor 

^junction = R+2r (5) 

and the leads are connected by an additional resistor 

Ri,^ = {R/r){R + 2r). (6) 

This is a consequence of the "Y-A" transformation of re- 
sistor networks"'^. For simplicity, we consider a capacitive 
coupling between the leads and central grain only and no 
Josephson coupling between the leads. Such a coupling 
will be generated in the renormalization flow. 



The imaginary-time effective action of the symmetric 
two-junction system can be written as a functional of the 
phase differences </>! and 02 across the first and second 
junction. 



Sc[(l>l,4>2\ + Sj[(t>i,(t>2] + SD[(t>l,4>2\-,{'^) 



where the charging term S'c, the Josephson coupling term 
Sj and the dissipation term Sd read 



'S'c[<^i,| 



WE^Jo ^"iidTi 

2] = -Ej / rfr[cos(0i) cos(02)], 
Jo 



(2) 
(3) 



Rq 



/ 

Jo 



drdr'- 



R{R + 2r)Jo sin2((7r//3)(r-T')) 
X [R{Mt) - Mr')? + R{Mr) - Mr')r 
+r{{Mr) + Mr)) - (Mr) + Mr')))'] ■ 

(4) 

Ec = e^/2C is the (single electron) charging energy of 
each junction and sets the overall energy scale. Ej de- 
notes the coupling strengths of the junctions. Ohmic 
dissipation in the resistors is introduced using the model 
of Caldeira and Leggett^^'^^. 



B. Weak coupling limit 

When the Josephson coupling energy Ej is small, it 
can be used as a small parameter for a perturbative RG 
analysis. As explained in Ref. 1, in addition to the bare 
Josephson energy of the two junctions, the RG flow pro- 
duces yet another Josephson coupling - the cotunneling 
J+. In the same sense that Ej is the amplitude for a 
pair-hopping between the leads and the grain, J+ is the 
amplitude for a Cooper-pair to tunnel between the two 
leads, skipping over the grain (see Fig. 3). To be more 
speciflc, the Cooper-pair hopping conductivities are 



Gab = Gbc ^ -E'j, Gac ^ J^- 



(7) 



Let us quote here the RG equations for the Josephson 
strengths in the symmetric case^ . To distinguish between 
the bare Jospehson energy Ej and the renormalized one, 
we use J to denote the flow of Ej. For the Josephson 
strengths J of the junction and J+ between the leads we 
get 



^ = J 1 



dJ 
H 
dJ+ 

~dr 



= J+ I 



R + r 

Rq 
2R 

Rq 



R TT 



Rc 



-J' 



(8) 
(9) 
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FIG. 4: Depending on the values of the resistors r and R, the 
low energy behavior of the two-junction system is described 
by three stable fixed points. In terms of the weak coupling 
Josephson variables J and J+, these points are: (i) Normal 
(NOR) fixed point at J = J+ = 0, where Josephson tunnel- 
ing is suppressed, (ii) Fully Superconducting (FSC) point in 
which both J and J+ are large and the flow is to strong cou- 
pling, (iii) SC* point where J+ flows to strong coupling and 
J = 0. Here, Cooper-pair tunneling is suppressed from lead 
to grain, but is very strong between the two leads. 



The Josephson coupHng J+ between the leads is origi- 
nally zero, but will be generated under the flow in the 
second order in J. From these RG equations we can infer 
the scaling behavior of the Cooper-pair conductivities in 
the asymptotic low-temperature regime and in the ICFP 
area. This is done by identifying the temperature T with 
the RG scale as follows: 



(10) 



The normal phase of the system is described by the 
fixed point J = J+ = 0, in which the Josephson junc- 
tions are insulating (sec Fig. 4). The Josephson coupling 
in this phase (and by Eq. (7) thus also the Cooper-pair 
conductivities) are expected to vanish as a power law in 
T, 



J. 



T 



(11) 



The signature of the normal phase is a drop of both lead- 
to-lead and lead-to-island conductance as temperature is 
reduced. 

If J+ is relevant at first order {R < Rq/2), then at low 
temperatures points A and C in Fig. 1 become short cir- 
cuited, so the leads become phase coherent independently 
of J. The RG equation for J becomes 



dJ 



J 1 



R/2 



Rq 



(12) 



When J is irrelevant (r -I- R/2 > Rq), the junctions are 
insulating despite the phase coherence between the leads. 



This is the SC*phase, and it is described by the fixed 
point shown in Fig. 4. The conductivities of the indi- 
vidual junctions in this phase vanish as the square of J, 
which follows the power law 



r + R/2 



(13) 



The lead-to-lead conductance diverges as a power law of 
the temperature and this divergence will be discussed in 
the next section using the strong coupling analysis. Note 
that the signature of the SC*phase is an increase in the 
lead-to-lead conductance, and a drop in the lead-to- island 
conductance as the temperature decreases. 

The FSC phase occurs when both J and J+ are rel- 
evant {r + R/2 < Rq). The weak coupling equations 
describe the initial increase in the Josephson couplings. 
The divergence of the conductance, however, can only be 
described in the next section, where the strong coupling 
limit is discussed. 



C. Strong coupling limit 

When Ej is large, one can no longer treat the Joseph- 
son coupling as a perturbation. But by employing a dual- 
ity, one can describe the system in terms of phase slips^ . 
Phase slips describe a sudden winding of a phase differ- 
ence across a junction. There are two kinds of phase slips 
we need to consider: individual phase slips and phase-slip 
dipoles. Individual phase slips create a potential drop, 
and hence dissipation in a junction independently of the 
other junction; we denote their fugacity by C- A phase 
slip in junction AB is partially screened if an anti phase- 
slip (a slip with the opposite winding) simultaneously oc- 
curs in junction BC. This partial screening leads to the 
generation of phase-slip dipoles in the RG flow, in a very 
similar fashion to the generation of J+. We denote the 
dipole fugacity by C^jj. By determining how often phase- 
slips occur, C and determine the effective resistance 
across the Josephson junctions (see Fig. 5): 

Rab = Rbc -C^ + Cl, Rac - 2C'. (14) 

The phase slip dipoles do not affect the lead-to-lead re- 
sistance since a dipole produces a voltage blip in one 
junction and an opposite voltage blip in the other junc- 
tion. Hence, the two blips add up to zero, and do not 
produce a voltage drop between the leads. 
The flow equations for C and are: 



dC 
dl 

dCo 
dl 



C(i- 

Cd ( 1 



RQ{R + r) 
2Rr + i?2 

2RRq 



RRq 
i?2 + 2Rr 

tRq ^2 



i?2 + 2Rr 7 i?2 -(- 2Rr 



CCd, (15) 
C- (16) 



In the FSC phase both ( and C,d are irrelevant, and 
decay to zero as a power law in temperature (and so do 
the resistances by Eq. (14)), 



Cd-T 



(17) 
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FIG. 5: In the strong-coupling limit we consider the RG flow 
of the phase slip fugacities ( and C,d ■ (a) ^ is the amplitude, or 
fugacity, for a single phase-slip across either of the junctions 
that connect the leads to the middle grain. When a phase 
slip occurs it produces a momentary potential drop across 
the junction, which leads to dissipation, (b) When r > 0, 
interactions between phase slips in the two junctions produce 
a phase-slip dipole in second order in the RG. The fugacity 
of the dipoles is C^d- When a dipole crosses the system there 
are two simultaneous momentary voltage drops with opposite 
directions on the two junctions, such that there is no potential 
drop between the leads, A and C. 



The signature of the superconducting phase is thus a 
drop in both the lead-to-lead and lead-to-grain resis- 
tance. This phase is controlled by the C = C-D — fixed 
point (Fig. 6). Note that the FSC and NOR phases are 
very similar in their behavior to the single junction prob- 
lem which has been investigated numerically in Rcf. 14. 

When (^D becomes relevant {r + R/2 > Rq), phase co- 
herence is destroyed between the leads and the central 
grain, but not between the two leads. At low tempera- 
tures, this practically implies that the effective resistance 
of the resistor r is diverging. Therefore the flow equation 
for ( becomes 



dl 



CI 



Rq 
2R 



(18) 



When R < Rq/2, even though (d has diverged, ( re- 
mains irrelevant and decays to zero at low temperatures 

as 



(19) 



In this phase the lead-to- lead resistance falls off as C^, 
but the lead-to-grain resistance diverges as 1/J^ from 
Eq. (13). Again, this is the SC*phase and the correspon- 
ding fixed point is shown in Fig. 6. 

In the normal phase, both C and Cd diverge at low 
energy scales, and one should use Eqs. (9) and (11) to 
describe the low temperature behavior of the measured 
resistances. 




FIG. 6: In terms of the strong coupling phase-slip fugacities 
and ("o, the three stable fixed points are the following: (i) 
Fully superconducting (FSC) fixed point sd, ( = = 0. In 
this point both single phase slips and phase-slip dipoles are 
suppressed, (ii) Normal (NOR) point in which both ^ and 
are large. There is no phase coherence and Cooper-pair 
tunneling in the system, which corresponds to weak coupling, 
(iii) The SC*point where (d fiows to strong coupling and ^ = 
0. Phase coherence is maintained between the leads (points A 
and C) but there is no phase coherence between the leads and 
the grain B. A quick look at Fig. 4 reveals the strong- weak 
duality of this system. 



D. Intermediate coupling fixed point (ICFP) 
region 

In the ICFP region a fourth fixed point appears (in 
addition to the FSC, SC*, and NOR fixed points), and 
the scaling behavior of J and J+, or ^ and Cd, is deter- 
mined at intermediate temperatures by this fixed point. 
At lower temperatures the asymptotic T scaling be- 
havior is determined by the J = J+ = (Normal phase) 
fixed point, or by the C = Cd = one (FSC phase). As 
opposed to the other regions of the phase diagram, where 
the resistors alone determine the phase, in the ICFP re- 
gion, the phase of the system also depends on the value 
of Ej / Ec ■ These parameters determine the critical value 
of the initial amplitude of Cooper-pair hopping, or, using 
the strong coupling picture, to the phase-slip fugacities. 
In this region there is no SC*phase, and the transition 
between the FSC and NOR phases is direct. An illus- 
tration of this situation in the ICFP region is given in 
Fig. 7. 

The existence of the unstable fixed point is evident 
from the nonlinear flow equations for the weak and strong 
coupling limits. We start with the weak coupling flow 
equations (Eqs. (8) and (9)) written as follows 



dJ 

H 

dJ+ 

~dr 



R 

~Ju + —JJ+, 
Uq 



(20) 
(21) 
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where 



R + r 
Rq 



1; W = — — 



2R 

Rq 



(22) 



When both u > and w > 0, the RG equations (20) and 
(21) have a third fixed point (in addition to zero and oo). 
This point is at 



J* = 



(23) 



Therefore the lines w, = and w = mark the weak- 
couphng boundaries of the ICFP region (see Fig. 2). 

Similarly, in the strong coiipUng limit we can write the 
flow equations (15) and (16) as 



dC 
dl 
dCp 
dl 



-Cu- 



RRr 



-TT = -Cdv + 



i?2 + 2Rr 

"^Rq /■2 



i?2 + 2Rr 



CCd, 



where 



RQjR + r) 
2Rr + i?2 



1, 



w 



2RRq 
i?2 + 2Rr 



- 1. 



(24) 
(25) 

(26) 



As in the weak coupling limit, when m > and W > 0, an 
unstable fixed point appears at intermediate values of Q 
and C_D • The lines u > and W > mark the boundaries 
of the ICFP region on the strong coupling side. The fixed 
point fugacities of the ICFP are given by 




i?2 + 2rR_ 
RqR 



-u. 



(27) 



The scaling properties of the system near criticality and 
at intermediate energy scales are determined by the crit- 
ical properties of the above ICFP and, in particular, the 
relevant and irrelevant directions and the exponents as- 
sociated with them: and A_. These properties are 
calculated in Appendix A. 

Qualitatively the situation is illustrated in Fig. 7. 
Consider a generic flow of C a-nd C/j at some value of 
the resistors r and R. If the system we are considering 
has initial fugacities and in the FSC-region, then 
both C and flow towards eventually, and the system 
is in the superconducting phase at low energies. But this 
behavior is not at all trivial: before decaying to zero, (d 
starts off growing as C decreases. Similarly, if and 
Cj)^ ~ are in the NOR region, once the flow passes by 
the ICFP, both and become relevant. But before 
the flow reaches the fixed point, there is an energy range 
in which (jj grows, but C decreases. Generally speak- 
ing, the extent of this range of energies is determined 
by the irrelevant critical exponent, A_ of the ICFP (see 
Appendix A). 

The FSC-NOR critical point occurs when C^°^ and c'^^ 
are on the critical manifold which flows exactly into the 
ICFP. In the two junction system = 0, and therefore 




FIG. 7: Typical flow of the parameters and ("d. If the initial 
values of and ("d arc in the region marked FSC they flow 
towards = ("u = after flowing towards the intermediate 
coupling fixed point. If they are initially in the region marked 
NOR, they flow towards the strong coupling and the FSC 
phase, again, after flowing to the ICFP flrst. The dashed lino 
is the critical manifold of the NOR-FSC transition, and the 
arrows near the fixed points mark the eigcndircctions of the 
linearized flow. The physical line is Co = 0, and therefore 
the critical Cc for the NOR-FSC transition is given by the 
intersection of the critical manifold with the Cd = axis. 



there is a critical fugacity above which the system is in 
the NOR phase, and below which it is in the FSC phase. 
The initial phase slip fugacity is determined by the ratio 
Ej/Ec, and therefore this transition (for a given value 
of the resistances R and r in the ICFP region) can be 
tuned by changing the Josephson energy. The qualitative 
picture described above is equally valid in weak coupling, 
where the only difference would be discussing the pair- 
tunneling amplitudes, J and J+, rather than the phase- 
slip fugacities. 



The flow of the fugacities (or pair-tunnel amplitudes) 

before reaching the ICFP determines the behavior of the 
resistance as a function of temperature at intermediate 
temperatures. In these temperature ranges the behav- 
ior of the resistance may be misinterpreted as any of the 
three phases of the system. Particularly, in the region of 
parameter space where ( decreases and Cd increases, the 
lead-to-lead resistance decreases, since it is only propor- 
tional to C^, but the lead-to-grain resistance increases. 
This behavior could be misinterpreted as the system be- 
ing in the SC*phase. In order to determine the true T = 
phase one has to investigate the system at very low tem- 
peratures. These crossover effects indeed appear explic- 
itly in the Monte Carlo simulations of the system. 
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III. MONTE CARLO METHOD 

Monte Carlo simulations of the two-junction model as 

shown in the lower part of Fig. 1 can be performed us- 
ing variants of the local and cluster updates detailed 
in Ref. 14. Imaginary time is divided into N time 
slices of size At = and the variables (/)i^„ and 4>2.n, 
n = 1, . . . , N (phase differences across the junctions at 
the discrete times t„ = nAr) are used to represent the 
phase configuration. We implemented the following types 
of Monte Carlo updates: 

1. Single junction cluster updates 

A cluster of connected sites is constructed in one 
of the junctions as outlined in Ref. 14. The cost 
in action of flipping the cluster, ASiead) which is 
associated with the last term in Eq. (4), 

~ i^icadio sin2((7r//?)(r-r')) 

(28) 

must be calculated and the cluster move accepted 
with probability 

p = min(l, exp(-A5iead))- (29) 

2. Single junction local updates 

Local updates in Fourier-space are proposed in one 
of the junction, as detailed in Ref. 14. They, too, 
are accepted with probability 

p = min(l, exp(-AS'iead))- (30) 

3. Compensated single junction cluster updates 

The SC*-phase is characterized by phase coherence 
between the leads (<?!'i(r) + 02 (t) « const) but 
strong fluctuations of the variables 4>i (t) and 02 (t) 
(insulating junctions). Therefore, the fluctuations 
in the two junctions essentially compensate each 
other and efficient updates in the SC*-phase should 
take this constraint into account. 

In a compensated single junction cluster update, a 

cluster of connected sites is constructed in one of 
the junctions as outlined in Ref. 14 and an equal, 
but opposite move proposed for each phase vari- 
able in the other junction. Since A0i -|- A(/)2 = 0, 
there is no cost in action associated with the 5icad- 
term in Eq. (4). However, there is a cost in ac- 
tion, A^oj, associated with the contribution from 
the other junction (hence the notation "oj"), 

-Rjunction Jo siu^ ((tt//?) (t - t')) 

(31) 



junction 1 jimction 2 




FIG. 8: Illustration of the two-junction cluster updates. The 

axis is chosen in at a symmetry point of the potential neigh- 
boring <t>i,n^^i^ (^^axis IS a random site), but on opposite sides in 
the junction i — 1 and 2. A second random site nroot is picked 
as the root site of the cluster. The sites connected to the root 
site are determined using the bond probabilities (33) and are 
marked with black dots. The now configuration is obtained 
by flipping the cluster around the axis. In other words, the 
spin like variables defined relative to the axis and shown as 
an arrow on the first site are inverted during a cluster update. 

The compensated cluster move should therefore be 
accepted with probability 

p = min(l,exp(-A5oj)). (32) 

4. Two-junction cluster updates 

A random site k is picked and an axis rii in each of 
the two junctions chosen among the two closest to 
4'i{Tk), such that niTT < and n27r > 02, fc or vice 
versa. Relative coordinates 4>f^^^ = (j^i — nf^^^n are 
introduced in both junctions and a cluster of sites 
connected to the root-site k is constructed using 
the bond probabilities 

p{k,l) = max(0, 1 - exp(-AS'fc,0) (33) 
where the cost in action of breaking a bond, ASk,i, 
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is defined as 



j=l,2 



axis 



xaxis 



)) 



;=i,2 



(34) 



In the above expression, g{j) is the kernel (j 7^ 0) 
1 



1 



Rq 



<5j,jv-i) 
(7r/iV)2 



8772 injunction sin((7r/iV)j)2 



(35) 



Hence, in contrast to the single junction case, a 
cluster can contain relative phase variables of both 
signs. The cluster building process (34) takes into 
account the capacitive, dissipative and Josephson 
contributions from both junctions, but not the dis- 
sipative contribution from the S'load-term. A two- 
junction cluster move therefore can only be ac- 
cepted with probability 

p = min(l, exp(-A5ioad))- (36) 



IV. PHASE DIAGRAM 

We first use the efficient Monte Carlo scheme outlined 
in Sec. Ill to identify the three phases NOR, FSC and 
SC*, and to determine the phase diagram for an inter- 
mediate value of the Josephson coupling Ej . This allows 
to test the theoretical predictions outlined in Sec. II and 
in Fig. 2. 

The result of this study is shown in Fig. 11 and ex- 
plained in Sec. IV A. The agreement between the Monte 
Carlo calculation and the theory is very good. Small 
deviations, however, appear in the vicinity of the (tri- 
critical) meeting point of the three phases: r = 0.75Rq 
and R = Rq/2. These deviations from the theoretically 
predicted phase diagram are explained in Sec. IV B as 
crossover effects, and are indirect evidence for the exis- 
tence of the intermediate coupling fixed point. 



A. Simulation results 

We use the resistance at imaginary frequencies to iden- 
tify the state of conductance between the leads and from 
the leads to the central grain. The imaginary frequency 
resistance is defined as 



Rq ^tt 



(37) 



where ujn = (27rn) / /3 denotes a Matsubara frequency and 
{4'4>)lj„ the Fourier transform of the phase-phase correla- 
tion function ((/)(0)(/)(r)) (0 = 0, in the case of conduc- 
tance from lead to central grain and <f> = <f>i + (j)2 for the 
conductance from lead to lead). 



r/fio=0.5 
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FIG. 9; Imaginary frequency resistances for r = 0.5 and 
l3Ec = 25, 50, 100, 200, 400, 800 and 1600. Only the first 
10 Matsubara points are shown for each temperature and the 
extrapolation of these curves to zero frequency gives the zero- 
bias resistance. From bottom to top, the different sets of 
curves correspond to i? = 0.5, 0.6, 0.625, 0.65, 0.675, 0.7 and 
0.75. The curves in the upper panel show the resistance from 
lead to central grain, those in the lower panel the resistance 
from lead to lead. In both cases the superconductor-to-normal 
phase transition occurs &t 0.65. 



In Figs. 9 and 10 we show the results of such an analysis 
obtained for Ej/Ec = 1 and I^tEc = 0.25. The inter- 
pretation of the data is the same as in the single-junction 
case discussed in Ref. 14. Figure 9 shows the transition 
across the FSC-NOR phase boundary at r = 0.5i?Q. The 
different sets of curves correspond to different values of 
R and the different lines in each set of curves to different 
temperatures. Only the lowest ten Matsubara frequen- 
cies are shown. Extrapolating these curves to w„ — > 
yields the zero bias resistance, which decreases to zero 
with decreasing temperature in the T = superconduct- 
ing state. If the junction turns insulating with decreasing 
temperature, the resistance increases and eventually sat- 
urates at the value r + R (from lead to central grain) or 
2R (from lead to lead). It is obvious from the data in 
Fig. 9, that the transition occurs simultaneously in the 
individual junctions and between the leads, in accordance 
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FIG. 10: Imaginary frequency resistances for r = 1.1 and 
PEc = 25, 50, 100, 200, 400, 800 and 1600. Only the first 
10 Matsubara points are shown for each temperature and the 
extrapolation of these curves to zero frequency gives the zero- 
bias resistance. From bottom to top, the different sets of 
curves correspond to R — 0.4, 0.45, 0.5, 0.55 and 0.6. The 
curves in the upper panel show the resistance from lead to 
central grain, those in the lower panel the resistance from 
lead to lead. For R < 0.5, the junctions are insulating, but 
the device from lead to lead is superconducting (SC* phase). 



with theoretical predictions. 

Figure 10 shows the transition across the SC*-NOR 
phase boundary at r=l.l. First of all, we note that the 
data for R < 0.5Rq clearly prove the existence of the SC* 
phase. While the individual junctions turn insulating as 
T ^ 0, the resistance from lead to lead decreases to zero 
as the temperature decreases (T = superconductivity). 
If R is increased, however, the device undergoes a (lead to 
lead) superconductor-to-metal transition at i? = O.bRq, 
as predicted by theory. 

Repeating this type of analysis for several values of r 
we could map out the phase diagram which is shown for 
Ej/Ec = 1 in Fig. 11. As expected, the FSC-NOR phase 
boundary lies somewhere in between the limiting values 
calculated analytically for Ej/Eq ^ 1 and Ej jEc <C 1, 
indicated by the dotted lines (see also Fig. 2). A close- 
up of the FSC-NOR phase boundary is shown in Fig. 12; 
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FIG. 11: Phase diagram obtained from the analysis of the 
resistance-versus-temperature behavior. We find a good 
agreement with theoretical predictions (black dotted lines), 
except in the region around the tricritical point, where the 
three phases meet. 
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FIG. 12: A closeup of the phase diagram in the ICFP re- 
gion for Ej/Ec ~ 1. The dashed line is the RG prediction 
for the FSC-NOR phase boundary as obtained in App. B, 
with the bare phase slip fugacity being the only fit parame- 
ter, determined using the critical point at r = Rq/^ to be 
^0 = 0.295. The RG-results agree well with the Monte Carlo 
results. Filled dots mark the FSC-NOR phase boundary and 
empty dots a FSC-SC*or SC*-FSC transition. 



the measured critical line at Ej / Ec = 1 agrees very well 
with an RG-based calculation (presented in App. B). 

Overall, we find a good agreement with the theoreti- 
cally predicted phase diagram, except in the vicinity of 
the tricritical point, r = 0.75Rq, R = Rq/2, where the 
three phases meet. Phase boundaries there appear to be 
shifted to a somewhat larger value of R and a smaller 
value of r. As it turns out, in this region, crossover 
effects are important and determining the phase from 
finite-temperature simulations can be misleading. A de- 
tailed discussion of the observed deviations is given be- 
low, in Sec. IV B. 

As the Josephson coupling strength Ej is increased 
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FIG. 13: Phase boundary as a function of the Josephson cou- 
pling strength calculated at v/Rq = 0.5. The dashed lines 
indicate the analytically predicted values for the critical R in 
the weak- and strong-coupling limits. 

from zero to infinity, tlie FSC-NOR pliase boundary 
sliifts from the weak coupling to the strong coupling limit 
indicated in Fig. 2. This behavior is demonstrated in 
Fig. 13, which shows a cut across the critical surface at 
t/Rq = 0.5. 

B. Crossover effects in the Monte Carlo results 

The theory of the two-junction system as outlined in 
Sec. II and in App. A leads us to expect crossover behav- 
ior at intermediate temperatures. Although the range of 
temperatures at our disposal is limited, and the crossover 
regime spans a narrow range of energy scales, we see in- 
direct evidence of the crossover effects in the measured 
phase diagram. As mentioned above, near the meeting 
point of the three phases, the predicted phase boundaries 
do not completely agree with the Monte Carlo simulation 
results (Fig. 11). We will explain these deviations using 
the RG analysis. 

The deviations of the measured data from the theo- 
retical predictions can be understood by looking at the 
RG-trajectories given by Eqs. (15) and (16). As ex- 
plained in Sec. II D, the flow of the phase slip fugacities 
or pair-tunneling amplitudes near the intermediate cou- 
pling fixed point can make the measured resistance as 
a function of temperature behave as though the system 
is in the SC*phase. Near the meeting point of the three 
phases, the flow towards the ICFP is extremely slow; it is 
dominated by the critical exponent A_ given by Eqs. (A9) 
and (A22): 

^(weak) ^ ^(-W- + 8uw) , (38) 

^(_strong) ^ 1 ^^__^_2_^_^^-^ ^ (39) 

where u, w and u, w are defined in Eqs. (22) and (26). 



Loosely speaking, but more specifically, in this regime 
X — X* ^ T^- with X being (, or Jj J+i and x* 
the respective fixed-point value. As can be seen from 
Eqs. (38), (39) and Fig. 2, A_ in both the weak and 
strong coupling regimes vanishes at the meeting point 
{w = W = 0), and therefore is expected to be small near 
this special point for any value of the Josephson cou- 
pling. This implies that since the Monte Carlo calcula- 
tion is limited to temperatures above T ~ Ec/2500 (if we 
choose AtEc = 0.25, which seems appropriate) it may 
not probe the ground state of the system, but rather the 
crossover physics. This will lead to distorted NOR-FSC 
phase boundaries in this region. As it turns out, this 
phenomenon is not restricted to the ICFP region: slow 
crossovers occur all around the triple point and may shift 
the observed NOR-SC*and SC*-FSC phase boundaries, 
as indeed is seen in Fig. 11. 

To demonstrate the above behavior we plotted the pre- 
dicted RG flow of the phase slip fugacities in the strong 
coupling limit for four values of r / Rq in the region of in- 
terest (Figs. 14 and 15). We color coded the plot accord- 
ing to the RG flow parameter Z = InA/T where A is the 
ultra-violet cutoff. Roughly speaking, the maximum RG 
scale we probe in the Monte Carlo calculation is ^max ~ 7, 
and therefore we stop the color coding at the value I = 10. 
The various flow lines all start with the same initial fu- 
gacity (same Ej/Eq), but have the resistance R varied 
across the transition. In the Monte Carlo simulation one 
measures the lead-to-lead resistance, which we expect to 
behave as 

Rac-C^ (40) 
and the lead-to-grain resistance, which goes as 

Rab = Rbc C + Cl- (41) 

As can be seen from Fig. 14a, at large values of r, 
r/Rg = 1.2, the flow within the accessible range of I 
allows an easy determination of the phase: at low ener- 
gies, C slowly decays or diverges to the left or the right 
of the black critical flow line at ii = Rq/2, while Cd 
strongly diverges (for this illustration we assumed that 
when (d > 0.7 the grain starts becoming effectively insu- 
lating). In Fig. 14b, at v/Rq = 0.8 closer to the meeting 
point (which is at r/Rq = 0.75) we see that this strong 
distinction cannot be made. There is a whole region of 
parameters to the right of the black line (marking the 
theoretically predicted critical flow), in which the lead- 
to-lead resistance {Rac ^ C^) decreases, and the single 
junction resistance {Rab ^ + Cd) increases. This re- 
gion eventually flows to the Normal fixed point, but in 
the finite temperature Monte Carlo simulations this can- 
not be observed and it appears that an SC* phase exists 
at values R > Rq/2, as shown in Fig. 11. To the left 
of the black line in Fig. 14b we theoretically expect the 
SC*phase, but even there we see a crossover which will be 
misinterpreted as an FSC phase: Rab decreases and, for 
a range of parameters, also Rac- Rac only starts to grow 
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FIG. 14: Simulated RG flows of phase slip fugacities 
(Eqs. (15) and (16)). The initial conditions for all plots 

are C*"' = 0.25, Cd' = 0. The color-coding marks the RG 
flow parameter I = InA/T, and regions not accessible with 
Monte Carlo simulations {I > 10) are shown in light gray, (a) 
r/Rq = 1.2, R/Rq varies from 0.4 to 0.6 from left to right. 
The black curve is at the critical value R — Rq/2. Within 
the energy scales probed by the Monte Carlo simulation it 
is easy to distinguish between the SC*and NOR phases. In 
the SC*phase C decreases making the lead-to-lead resistance 
Rac ~ decrease at low energies, while grows rapidly, 
producing an insulating regime for the effective resistance for 
each junction, Rab- We assume that begins to saturate 
at (d = 0.7 and makes the flow cross over to Eq. (18). (b) 
r/Rq = 0.8, and R/Rq varies as in (a). At this value of r 
the temperature range accessible with simulations is no longer 
sufficient to determine the true phase boundaries. For a range 
of 7? values near and above the critical point R = Rq/2 we 
may mistakenly identify the phase as SC*, since Rac ~ 
is decreasing, but Rab ~ + Cd is increasing. Similarly, 
for R < Rq/2 it looks as if the system is in the FSC phase 
(instead of the SC*phase), since (d grows very slowly. This 
explains the deviations of the numerically determined phase 
boundaries (Fig. 11) from the theoretical ones. 



at a lower energy scale, disclosing the true SC*phase. In- 
deed, in Fig. 11 we see that the Monte Carlo calculation 
indicates an FSC region at t/Rq — 0.8 and R < 0.5Rq, 
where it should be the SC*phase. 

In the ICFP region we again encounter slow crossovers 
associated with the intermediate coupling fixed point. As 
can be seen in Fig. 15c with r/Rg = 0.7, again there 



r/Rq =0.7 
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FIG. 15: Same as Fig. 14: (c) r/Rq = 0.7, R/Rq varies 
between 0.5 and 0.6 from right to left. These flows are in the 
ICFP region. Even though in this region we expect a direct 
NOR-FSC transition, at the energy range accessible in the 
Monte Carlo simulation we again see a region which would be 
mistaken as SC*, which is indeed the case in Fig. 11. The 
slow crossover is due to the proximity to the triple point at 
r — 0.75Rq where the three phases meet, (d) r/Rq — 0.5, 
and R/Rq varies as in (c). As in Fig. 14a even within the 
energy range accessible numerically, we clearly see the NOR 
and FSC phases showing up - and (d either both diverge or 
both decrease at the lowest temperatures accessible by Monte 
Carlo simulation. 



is a range of parameters which would be mistaken for 
the SC*phase, in which Rac seems to drop, while Rab 
grows. Since the RG is stopped approximately at the 
^max corresponding to our Monte Carlo calculation, the 
observed flow in this case is dominated only by the critical 
exponent A_ of Eqs. (38) and (39), which vanishes at 
the triple point. For comparison, in Fig. 15d we show 
the RG flow for r/Rq = 0.5. In this case the RG does 
flow to the stable NOR and FSC fixed points at energy 
scales higher than the lowest temperature accessible in 
our calculation. Indeed, in this parameter range there is 
no longer any evidence of an SC*-like phase in the Monte 
Carlo calculation. 
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V. THE CRITICAL TWO-JUNCTION SYSTEM 
IN THE ICFP REGION - COMPARISON WITH A 
SINGLE JOSEPHSON JUNCTION 

In this section we investigate the direct NOR-FSC 
transition, and compare the critical behavior of the two- 
junction system with that of a single junction with the 
same Ej/Ec- We make a rather surprising observation: 
along this phase boundary, the effective resistance of the 
junction, the temperature dependence of the mean phase 
fluctuations, and the correlation exponents are within 
error-bars the same as those in a single resistively shunted 
Josephson junction at criticality (with the same Ej/Ec 
and AtEc)- This suggests that many of the features 
of the NOR-FSC transition at the ICFP can be under- 
stood in terms of the single junction Schmid transition, 
although the ICFP is an interacting fixed point. We will 
first present the numerical results and then proceed to 
discuss them in Sees. VB and VI. 



A. Numerical results and resemblance to the single 
junction 

Associated with the drift of the phase boundary in 
the ICFP region as a function of Ej/Ec (see Fig. 13) 
is a continuous change in the critical resistance of the 
junctions and in the value of the correlation exponents 
(defined below in Eq. (42)). As mentioned above, these 
features of the fixed point are remarkably similar in the 
two-junction model and in a single resistively shunted 
junction at criticality. To illustrate this, we first consider 
the critical resistance and plot in Fig. 16 the resistance 
from lead to central grain as a function of inverse tem- 
perature for Ej/Ec = 1, r = Q.5Rq and R = 0.65Rq, 
which is a point on the FSC-NOR phase boundary. Also 
shown in the figure is the temperature dependence of the 
resistance in a single junction (with Ej/Eq = 1 and the 
same discretization step AtEc — 0.25) for several values 
of the shunt resistance Rs ■ The curves from top to bot- 
tom correspond to Rq/Rs = 0.9, 0.95, 0.975, 1.0, 1.025, 
1.05 and 1.1, respectively. For Rq/Rs > 1, the junction 
turns superconducting as T ^ 0, whereas for Rq/Rs < 1 
it becomes insulating. At the critical point Rq/Rs — 1, 
the resistance of the system (junction plus shunt resistor) 
is precisely the same as the critical resistance from lead 
to central grain in the two-junction model. 

By and large, the critical resistance of the two-junction 
system does not vary along the FSC-NOR phase boun- 
dary for fixed Josephson coupling. We plot its value in 
Fig. 17 as a function of r and compare it with the corres- 
ponding result for a single resistively shunted junction. 
A very good agreement is evident for t/Rq — 0.1, 0.2, 
and 0.5. For t/Rq = 0.7, closer to the meeting point 
of the three phases, the agreement is less good although 
still within error bars. 

Another quantity we measured is the mean phase fiuc- 
tuation ((0—0)^), which we found in Ref. 14 to grow pro- 
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FIG. 16: Critical lead-to-grain resistance as a function of in- 
verse temperature for Ej/Ec = 1 and AtEc = 0.25. The 
point r — 0.5Rq, R = 0.65Rq on the FSC-NOR phase boun- 
dary has been selected (see Fig. 11). For comparison, we also 
plot resistance- versus-temperature data for a single resistively 
shunted junction with Ej/Ec = 1 and AtEc = 0.25. The 
values Rs of the shunt resistors are (from top to bottom) 
Rq/Rs = 0.9, 0.95, 0.975, 1.0, 1.025, 1.05, and 1.1, respec- 
tively. At the phase transition point Rq/Rs = 1, the resis- 
tance of the single junction is exactly the same as the critical 
lead-to-grain resistance in the two-junction system. 
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FIG. 17: Value of the critical lead-to-grain resistance along 
the FSC-NOR phase boundary for Ej/Ec = 1. The critical 
resistance is nearly independent of r and within error-bars 
it is the same as in a single junction. The solid line shows 
the critical resistance of the single resistively shunted junction 
with Ej /Ec = 1 and Rs = Rq ■ Error estimates are indicated 
by the dotted lines. The dashed line is the result of an RG- 
based calculation of the effective lead-to-grain resistance. 



portional to the logarithm of the inverse temperature at 
criticality. The same is true in the two-junction system. 
In Fig. 18 we plot ((</)— c^)^) as a function of inverse tem- 
perature. The lines show the data for a single junction 
with the bold line marking the logarithmic growth of the 
phase fiuctuations at the critical point. The diamonds, 
circles and triangles show the data obtained from the two- 
junction system at criticality. They correspond to points 
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FIG. 18: Phase fluctuations ((0 — 0)^) as a function of inverse 
temperature. The curves from top to bottom show the data 
for a single resistively shunted junction with Rq /Rs = 0.9, 
0.95, 1.0, 1.05, 1.1, and 1.2 respectively. The symbols show 
the same quantity calculated in the two-junction model at 
criticality. Three points along the FSC-NOR phase boun- 
dary have been chosen and all the data agree with those of 
the single junction with identical values of Ej/Ec = 1 and 
ArSc = 0.25. 
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FIG. 19: Correlation function (42) for q = 0.25 and l3Ec = 
1600. The curves correspond to r = 0.5 Jig, and from top to 
bottom: Ej/Ec = 2, 1, 0.5 and 0.25, respectively; R in these 
plots is tuned to the FSC-NOR critical line. The solid lines 
show the exponents computed for the two-junction system 
and the dashed lines those obtained for a single resistively 
shunted Josephson junction at criticality. The dotted lines 
indicate the error on the two-junction calculation originating 
from the uncertainty on the critical resistances. 



along the FSC-NOR phase boundary with Ej/Ec = 1 
and AtEc = 0.25, as in the single junction case. The 
uncertainty on the critical value of R has not been taken 
into account in the error bars of Fig. 18. This uncer- 
tainty can account for the small deviations between the 
different data sets. 

In addition, we considered the exponents ri{q), which 
we measured for several values of Ej/Ec at r = Q.5Rq 
and with R/Rq corresponding to the FSC-NOR phase 
boundary. We define the latter as in Ref. 13 using the 
correlation function 

corr,(T) = {eMiQiHr) - <^(0))]), (42) 

with q some non-integral real number and (p the phase dif- 
ference across one of the junctions. In the normal phase 
and at criticality the correlations decay as 

corr,(T) ~ r-^''^'^^ (43) 

In the superconducting phase one observes a power-law 
decay of the connected correlation function 

corr,(r) = (exp[ig(</.(T)-0(O))])- |(cxpM(t)])|2. (44) 

In Fig. 19 we plot corrg^i/4(T) for the Josephson coupling 
strengths Ej/Ec = 0.25, 0.5, 1, 2 (as in Fig. 13) and 
the corresponding critical resistances r = 0.5 and R = 
0.5375(200), 0.5875(150), 0.650(20), 0.6875(200), respec- 
tively. We compare these correlation functions to those 
obtained for a single jimction at criticality (i?shunt = Rq) 
for the same values of Ej / Ec and the same discretization 
step AtEc = 0.25. As can be seen in Fig. 19, the corre- 
lation functions, and thus also the critical exponents r]{q) 
perfectly agree for Ej/Eq = 2 and 1, while they agree 



within error bars for Ej/Ec — 0.5 (the dotted lines in 
the figure show the correlation functions computed at 
R = i?criticai ± error). For Ej/Eq = 0.25 the agreement 
is no longer as good, but the exponents for a single junc- 
tion and for two junctions remain close and well within 
error-bars. 

While the exponent ri{q) varies as a function of Ej/Ec 
in a remarkably similar way as in a single junction, there 
is no dependence of this exponent for fixed Ej/Ec on 
the values of r and R, that is on the position along 
the FSC-NOR phase boundary shown in Fig. 11. Re- 
member that also the effective lead-to-grain resistance 
was within error-bars independent of r on the FSC-NOR 
phase boundary (Fig. 17). 

We plot 7?sin(9), obtained from a fit to 

^^"^ = sin((V^)r)^'^3.- ' ^''^ 

along the FSC-NOR phase boundary for = 0.5, 13 Ec = 
1600 and Ej/Ec = 1 in Fig. 20. The measured ex- 
ponents range between 0.05 and 0.06 and are thus well 
within error-bars. To estimate the error, which mainly 
originates from the uncertainty of ±0.02i?Q on the criti- 
cal value of Rc (for a given r), we also plot the exponents 
measured at Rc ± Q.Q2Rq in the figure. 

B. RG analysis of the critical effective resistance in 
the two-junction system 

The extensive investigation in the previous Sec. VA 
demonstrated a remarkable resemblance between the be- 
havior at criticality of the two-junction system and a sin- 
gle junction with the same Ej / Ec- In this section we will 
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FIG. 20: Correlation exponents r;sin(g = 0.5) obtained from 
a fit to /(r) = a/ smiTT/PTf""'" along the FSC-NOR phase 
boundary (see Fig. 11) for constant Ej/Ec ~ 1 and l3Ec ~ 
1600. The main source of error is the uncertainty of ±0.02_Rq 
on the critical resistance Rc for a given value of r. The in- 
dependence of the exponent rjsin{q — 0.5) with respect to r 
and the constant lead-to-grain resistance shown in Fig. 17 in- 
dicate that the critical properties only depend on the value of 
Ej/Ec. 



account for this resemblance by comparing the measured 
effective resistance at criticality of the two-junction sys- 
tem with the resistance as predicted from the RG flow 
equations (15) and (16). We emphasize that such an 
analysis is only approximate, since Eqs. (15) and (16) 
are derived in the limit of large Ej /Ec , and the case we 
are concentrating on is that of an intermediate Ej/Ec- 
Nevertheless, the strong-coupling RG equations provide a 
rather good fit to the behavior of the system at criticality. 
Since all quantities considered in Sec. V A (lead-to-grain 
resistance, phase fluctuations, and phase correlations) are 
related to the effective lead-to-grain resistance, the ana- 
lysis of the latter sufflces. 

The first step in the analysis is to connect the measured 
lead-to-grain effective resistance with the resistances R 
and r, and with the phase-slip fugacities ( and (d (which 
we will later substitute with their fixed-point values, C* 
and Ch)- This amounts to solving the circuit shown in 
Fig. 21. The junctions in the figure were each replaced 
by two components: a resistor 

Rc ~ aC\ (46) 

and a voltage drop due to phase-slip dipoles, Vd- The 
latter is given by 

Vo^aClih-h). (47) 

The solution of the circuit in Fig. 21 is sketched in 
App. C. It leads to the following expression for the lead- 
to-grain resistance 



R. 



eff ■ 



iR/a + C^){2r + R + aC 




FIG. 21: A circuit describing the measurement of the lead- 
to-grain resistance. Each junction is replaced by two compo- 
nents: a resistor 7?^ = aC,^ which describes the response of 
single phase-slips, and a voltage drop Vd which in junction 
BC is qC|)(-^i — I2) and in junction AB is ck^|j(/2 — h)- The 
Ohm-meter is considered as a current source, I. 



In Eqs. (46) and (47) we used the coefficient a to con- 
vert the square of the phase-slip fugacities to resistances. 
This coefficient is going to be the only input into the 
RG fiow analysis of the effective lead-to-grain resistance. 
An estimate for this coefficient can be obtained from the 
data point at Ej/Ec = 1 and r/Rg = 0.5 in Fig. 17, in 
conjunction with the measured single-junction effective 
resistance. The critical shunt resistance at this point is 
R/Rq = 0.65 ± 0.02. From the integrated RG fiow we 
can find Co such that for r = 0.5i?Q and R = 0.65i?Q the 
point (Co, 0) lies on the critical manifold. We find 



Co = 0.295 ±0.01. 



(49) 



This Co corresponds to Ej/Ec = 1- Note that Co is a 
function of Ej/Eq only in the strong coupling limit, and 
should only depend weakly on the resistances r, R (c.f. 
Rcf. 16). The same applies to the coefficient a. Bearing 
this in mind, we can calculate a by using the measured 
resistance of the single-junction system. This measured 
resistance is the resistance of the junction, aC^i paral- 
lel to the shunt resistor, which equals Rq at criticality. 
Thus, for Ej /Ec = 1 we have (Fig. 17) 



aCo ■ Rq 
Rq + aCo 



0.265 ±0.015 



from which we find 



a = 4.15 ±0.35. 



(50) 



(51) 



(48) 



Now we have all the pieces to predict the effective lead- 
to-grain resistance, Ros on the FSC-NOR critical line 
Ej/Ec — I- By using Eq. (48) with the fixed-point 
values C* and Czj from Eqs. (27), and a from Eq. (51), 
we obtain the curve of i?off as a function of r. The re- 
sult of this calculation is shown in Fig. 17 by the dashed 
line. One can see that the resistance predicted by the 
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FIG. 22: Illustration of the mean-field theory for the FSC- 
NOR transition in the two-junction model. A Josephson junc- 
tion at criticality is replaced by a resistor Rj, such that the 
resistance of the circuit on the top right corresponds to the 
measured resistance of the device. Replacing one of the junc- 
tions in the two-junction system by Rj, one can calculate the 
effective shunt resistance R^^ of the other junction. Identify- 
ing Rl^ with Rq yields an equation for the FSC-NOR phase 
boundary. 



RG changes little in the entire range, and remains in rea- 
sonably good agreement with the observed Monte Carlo 

ReS- 

Note that we assumed that the measured resistance is 
due to the fixed point characteristics of the two-junction 
system. Prom Sec. IV B, however, we know that in the 
vicinity of r = 0.75i?Q, where the three phases meet, 
crossover effects are dominant. Therefore we expect the 
lead-to-grain resistance calculated in this section to de- 
viate from the measured i?cfr in that vicinity. Another 
caveat for the current calculation is that it is correct 
up to second order in the phase-slip fugacities; fourth- 
order contributions to i?off arc neglected (although we 
keep fourth order terms in Eq. (48)). These corrections 
may also account for deviations from the measured Res- 



VI. THE ICFP AS A SELF-CONSISTENT 
FIXED POINT 

In the previous section we investigated the FSC- 
NOR transition extensively, and demonstrated a remark- 
able resemblance between the two-junction and single- 
junction systems at criticality. We adequately explained 
this surprising resemblance using the RG from Sec. II. 
In this section, however, we provide yet another explana- 
tion, albeit ad-hoc, for this resemblance. The alternative 
explanation is that the ICFP can be approximated as a 
self-consistent fixed point. In such a mean- field theory, 
illustrated in Fig. 22, the physics of a single junction 
emerges naturally. The idea behind this approach is that 
a Josephson junction undergoes a SC-NOR transition 
when its effective shunting resistance is Rq. 

Consider the two-junction system with phase slip fu- 
gacities Co (with = 0) in both junctions. In mean- 
field, one junction sees the other as an effective resistor 
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FIG. 23: Upper figure: SC-NOR phase boundary obtained 
from the mean-field prediction Eq. (55) for Ej/Ec ~ 1. The 
critical resistance Rc of the single junction with Rs = Rq 
is Rc/Rq w 0.265 (see Fig. 17). The dotted lines show how 
the error on Rc influences the outcome. A remarkably good 
agreement with the Monte Carlo results is obtained. 
Lower figure: A closcup of the phase diagram in the ICFP re- 
gion for Ej/Ec = 1 (Co = 0.295). The solid line is the mean- 
field phase boundary from Eq. (55) and the dashed line is the 
RG prediction for the FSC-NOR phase boundary. The RG- 
results agrees well with the mean-field prediction and with 
the Monte Carlo results for r < 0.6Rq. Filled dots show the 
position of the FSC-NOR transition, open dots indicate a SC- 
SC* or SC*-FSC transition, respectively. See section IVB for 
a discussion of crossover effects. 



(Figs. 22 and 24) with resistance 
Rj = aCo> 



(52) 



where Rj is defined from Eq. (50). Therefore the effective 
shunting resistance on each of the junctions is 



Ra ~ R 



r [R + Rj) 



(53) 



r + R + Rj 
and criticality is obtained when 

Rs = Rq- (54) 
From Eqs. (53) and (54) we find an expression for the 
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FSC-NOR phase boundary, 

^ {R + Rj)(Rq-R) 
^ 2R + Rj-Rq 



R 



(55) 



as a function of the parameter Rj, which depends on 
Ej/Ec- The limiting cases for strong and weak Joseph- 
son potential work out correctly. For Ej/Eq oo, 
Rj (superconducting junction) and Eq. (55) reduces 
to 



R{Rq - R) 



(56) 



whereas for Ej/Ec 
and we obtain 



2R-Rq 
0, Rj ^ oo (insulating junction) 



r = RQ-R. 



(57) 



These are indeed the equations describing the FSC-NOR 
phase boundaries shown in Fig. 2. 

In the mean-field approximation, each junction at crit- 
icality behaves as an independent junction, and thus ex- 
hibits the same phase correlations and fluctuations as the 
single junction with the same Rj, and hence exhibits the 
same dependence on Ej / Ec ■ 

Although this explanation seems naive, it does allow 
a simple and quantitative understanding of the observed 
(and calculated) properties of the two junction-system 
at the SC-NOR transition point. As further evidence, in 
Fig. 23, we compare the phase boundaries in the ICFP 
region as obtained by the simple mean-field treatment 
(Eq. (55)), the RG calculation (explained in App. B), 
and the Monte Carlo results. For r/Rq < 0.7 there is 
a good agreement between the three approaches. In the 
vicinity of the meeting point of the three phases, the RG 
boundary traces the measured FSC-SC*phase boundary, 
while the mean field prediction agrees with the Monte 
Carlo SC*-NOR phase boundary, even for values of r > 
0.75Rq. Note, though, that despite the agreement with 
the Monte Carlo results, the mean-field argument does 
not apply to the SC*-NOR phase boundary. 

The good agreement of the mean-field boundary with 
the measured FSC-NOR boundary seems rather surpris- 
ing. It can be explained once we understand why the 
mean-field treatment is (at least qualitatively) reason- 
able. The mean-field approach makes two assumptions: 
(1) a junction with phase-slip fugacity C can be replaced 
by a resistor a('^, and (2) a(1) can be neglected compared 
to r (see Fig. 24). At this point we can only justify the 
first assumption on qualitative, but not on analytical, 
grounds. Assumption (2) on the other hand is justified 
by the fact that the fixed point values of the phase-slip 
fugacities obey a{(^'^)^ <C r for much of the ICFP region. 
For Co = 0.295 (corresponding to Ej/Eq = 1) this ratio 
rises slowly from zero at r = to 0.25 at r = 0.6i?Q; it 
crosses 1 at r = 0.73i?Q. In the region where a{C,'^Y is 
no longer negligible, i.e., r > 0.6i?Q, crossover effects set 
in, and the RG flow is cut off at the lowest energy scale 
of the Monte Carlo simulation. In this case (^d does not 
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FIG. 24: Effective circuit in tiie mean-field approximation - 
taking phcise-slip dipoles into account. Junction AB is ap- 
proximated by a resistor with resistance = aC,"^ which is 
due to single phase slips. In addition, phase-slip dipoles add 
Rd = qCd to the resistor r. The location oiRo is determined 
such that it responds to the difference in the current on the 
two junctions (c.f. Eq. (47)). This circuit is equivalent to the 
one shown in Fig. 21 with J = 0. In the mean-field approxi- 
mation we assume that a NOR-FSC phase transition occurs 
when the effective shunting resistance that a phase slip on 
the junction BC probes equals Rq. Furthermore, we assume 
that Rd is negligible - an assumption which is well justified 
with respect to the Monte Carlo calculation in the ICFP re- 
gion and near the meeting point of the three phases, where 
crossover effects are evident (c.f. Fig. 22). 



reach its fixed point value and will be observed as small. 

The crossover effects justify assumption (2) of the mean- 
field in the same way that they explain the deviations 
from the predicted phase diagram in Sec. IV B. Thus 
the mean-field approach reproduces the phase diagram 
with cross-over effects. 

Above we made the assumption that in the mean-field 

approach a junction undergoes the SC-NOR transition 
when the effective shunting resistance sees equals Rq. 
We also derived the phase boundary using the RG fiow 
(App. B) without making any such assumption. This 
raises the question: how close to Rq is the effective shunt- 
ing resistancx; along the RG-predicted phase boundary? 
This effective shunting resistance is given by 



R 



leff , 



R + 



{R + Rc){r + RD) 
r + R + Rc + Ro ' 



(58) 



which can be deduced by observation from Fig. 24, with 
Rd = aCi) and i?^ = a('^. In Fig. 25 we plot this effective 
shunting resistance for Co = 0.295, which corresponds to 
Ej/Eq = 1, as a function of r for the entire ICFP re- 
gion, < t/Rq < 0.75. We assume that ( and (d take 
their fixed point values from Eqs. (27). As can be seen, 
the effective shunting resistance stays very close to Rq 
over the entire range. At this point we are unable to say 
whether this is just a coincidence (as the success of the 
mean- field treatment above may therefore be), or a uni- 
versal property of the phase boundary, which the approx- 
imate RG roughly reproduces. To answer this question, 
a more thorough consideration of the sine-Gordon model 
is necessary. 
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FIG. 25: The RG-based calculation of Rf (solid line) given 
in Eq. (58) as a function of r/Rg on the phase boundary in 
the IGFP region with Ej/Ec = 1 (Co = 0.295). The effec- 
tive shunt resistance is very close to Rq in the entire region, 
giving more validity to the proposed mean-field approach to 
the double sine-Gordon model of the two-junction system. 
The dotted line is the effective shunt resistance when we set 
R( — and Rd = 0, and is given for reference (this is just 
R + rR/{r + R)). 



VII. CORRELATION FUNCTION IN THE NOR 
PHASE 

In this section we compute and briefly consider the 
decay of the correlation function (see Eq. (42)) 

corr,(r) ee (exp[iq{(j){T) - 0(0))]) 

in the NOR phase of the system. We set r = 0.5Rq and 
R = 0.65Rq, and calculate corrg(T) for different values 
of the Josephson coupling energy. This choice of r and 
R corresponds to a critical value (-Ej/i?c )ciiticai — 1 (see 
Fig. 13). Our results are shown in Fig. 26. Interestingly, 
each curve has an approximate power-law behavior, but 
with exponents that vary with Ej/Ec- The correlation 
exponents 77sin(9 — 0.5) obtained from a fit to Eq. (45) 
are shown in Fig. 27 for the temperatures l3Ec = 1600, 
800 and 400. This behavior is very different from that 
seen in a single resistively shunted junction, where differ- 
ent Ej/Ec (at fixed shunt resistance Rs) yield the same 
exponents in the NOR phase^'*. 

However, one point should be kept in mind: even 
though in the previous section, we demonstrated and dis- 
cussed the similarities between the two-junction model 
and the single junction at criticality, there are also stark 
differences between the two systems. While in the sin- 
gle junction, flow lines corresponding to systems with 
the same shunt resistance but different Ej /Ec merge (as 
they flow parallel to the Ej/Ec axis only), in the two- 
junction system, under the same circumstances, they do 
not overlap (see Fig 7). Indeed, in Eqs. (9) and (16) we 
need an extra parameter, J_|_ and <^d, in order to write 
the renormalization group flow equations for the two- 
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FIG. 26: Correlation functions corrg=o.5(r) in the FSC phase 
for different temperatures and values of Ej/Ec, but fixed 
resistance r = 0.5Rq and R — 0.65Rq (corresponding to 
(i5j/_Bc')criticai = !)• Different Ej/Ec yield considerably dif- 
ferent exponents ??(?). 
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FIG. 27: Correlation exponents 77sin(g = 0.5) obtained from 
a fit to Eq. (45) for the temperatures f3Ec = 1600, 800 and 
400. 



junction system. 

In the single-junction model we expect systems with 
different J oc Ej /Ec but same shunt resistance to show 
the same low energy behavior. More specifically, at weak 
J we expect 



corrg (t, J) — b corr^ 



J 



(59) 



To compare two different systems, with junctions Ji and 
J2, we can choose 



i/(i?yflQ-i) 



(60) 
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and obtain 



R/Rn=0.65 



r/Rn=0.5 



corr, (r, Ji) = 



+ 



fJl 



l/(fl/flQ-l) 



(61) 



If the decay of the correlations is a power law, corrg (r) ~ 
|t-|-2i)(9)^ then the scaling relation (61) clearly shows that 
J must drop out of the power law. In fact, from Eq. (59) 
we find directly 



r]{q) = (f 



R 



(62) 



In the two-junction system, we observe a power law 
decay of corr^ (r, Ji), but with 77 (q) which seemingly de- 
pends also on J. In the absence of our understanding of 
the RG flows of the two-junction system, we might guess 
that a scaling relation similar to Eq. (59) holds, but with 
the shunt resistors, R and r also flowing. This possibility, 
however, seems to contradict Ref. 3, which claims that 
in sine-Gordon models the resistance can not get renor- 
malized due to its singular |a; [-frequency dependence. 

We can, however, qualitatively explain the dependence 
of rj{q) on J in terms of the RG flow equations (21). In 
terms of the two-junction RG flows for J and J+ (or C 
and Cd), a scaling form such as Eq. (61) does not hold 
in the two-junction system in the ICFP region. Instead, 
we need to integrate Eqs. (20) and (21). In Fig. 28, the 
RG flows corresponding to i? = 0.65i?Q, r — 0.5Rq, 
and Jo = 0.5, 0.4, 0.3, 0.2 {Ej/Ec = 1, 0.8, 0.6, 0.4), 
are shown. In the accessible energy range, all flows take 
different paths in the (J, J+) space. Also, they are all 
far from saturating at J = J+ = 0. The latter obser- 
vation provides the reason for the slower decay of the 
correlations for higher J: in the accessible energy range 
the two-junction system seems to have a larger effective 
J and J_|_, the higher the bare Jq is. These finite J and 
J+, in tmm, provide extra conductance in the system, 
and thus strengthen phase coherence, and suppress the 
decay of correlations. 

We furthermore note that the fiow in the temperature 
range 400 < PEc < 1600 (corresponding to 6 < / < 7.4) 
is rather slow, which might explain why there is hardly 
any temperature dependence in the measured rj (Fig. 27). 
However, we cannot rule out a scenario where the flows 
have actually converged to a set of distinct fixed points 
(with e.g. different values of J+), although the origin of 
such additional fixed points remains unclear. 



VIII. CONCLUSIONS 

In this paper, we presented an RG analysis and simu- 
lation results for a symmetric two-junction system. The 
phase diagram was found to contain all three phases ex- 
pected by theory: (i) fully superconducting (FSC), (ii) 
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FIG. 28: The integrated weak-coupling RG flows of Eqs. (21) 
in the (J, J+) plane for 7? — 0.65Rq and r — 0.5Rq. The 
value Jo = 0.5 corresponds to the critical value Ej/Ec = 1. 
These flows show qualitatively how the paths with different 
bare coupling Jo essentially do not overlap, and also do not 
reach zero in the accessible energy range. The flows are color 
coded to show the evolution of the logarithmic RG scale I = 
In Ao/A along the flow (Aq is the bare UV cutoff' and A is the 
renormalized cutoff). In Fig. 26 Z < 8, as in this plot. 



normal (NOR) and (iii) superconducting from lead to 
lead but normal from lead to grain (SC*). Furthermore, 
the phase boundaries obtained by the Monte Carlo sim- 
ulation agreed with the theoretical prediction, except 
for deviations near the tricritical point where the three 
phases meet. These deviations were interpreted as finite 
temperature crossover effects, and give indirect evidence 
for the existence of the intermediate coupling fixed point. 
The behavior of correlation functions in the NOR phase 
may also be indirect evidence for the flow of the coupling 
J+, as predicted by the RG-theory. 

The two-junction system possesses the unusual ICFP 
region, in which the NOR-FSC transition can be tuned 
by Ej/Ec as well as by the dissipation. We calculated 
this phase boundary using the RG treatment of Ref. 1, 
and found that it fits the Monte Carlo result quite well 
(c.f. Fig. 12). Theoretically it was shown that in this re- 
gion the properties of the two-junction system are deter- 
mined by a new intermediate coupling fixed point, with 
continuously varying exponents. 

An efficient Monte Carlo algorithm was used to tho- 
roughly investigate the intriguing ICFP region and the 
NOR-FSC phase transition. We calculated the lead-to- 
grain resistance, average phase fiuctuations, and phase 
correlations for several locations on the critical NOR- 
FSC manifold. This manifold is two dimensional, and 
can be parameterized by r and Ej / Eq ■ 

For the various sets of (r, Ej/Ec) on the NOR-FSC 
critical manifold we made a surprising observation: the 
effective fixed point lead-to-grain resistance, Rab, as well 
as other critical properties of the two-junction model, 
were nearly identical to those of a single junction (with 
the same Ej/Ec) at criticality {Rs = Rq). Moreover, 
the effective lead-to-grain resistance is within error-bars 
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independent of r and R, and depends only on Ej/Ec- 

Given the success of the RG analysis of the two- 
junction model, we tried to use it in order to account for 
the above observations. Using simple scaling arguments, 
we extended the RG of Ref. 1 and Sec. II to compute 
the effective lead-to-grain resistance at the critical point. 
This calculation also yielded a nearly r independent ef- 
fective resistance, and since in the limit r = 0, R = Rq, 
the two-junction system and the single junction system 
coincide, the observations about the resemblance of the 
two systems seem to be satisfactorily explained. The 
approximate RG-bascd resistance calculation fitted the 
QMC results quite weU (Fig. 17). 

But the observed similarity of the single- and two- 
junction systems could also be qualitatively interpreted 
in a more interesting way. The observation can be taken 
to indicate that at criticality, each of the two junctions in 
the system sees an effective environment that imitates an 
Rs = Rq shunt, and therefore shows the same effective 
fixed point resistance and critical properties as a single 
isolated junction at criticality. This mean-field argument 
was developed in Sec. VI, and it allowed us to predict 
the location of the FSC-NOR phase boundary at interme- 
diate Josephson energy well within error bars (Fig. 23). 
In addition, we used the RG predictions to compute the 
effective shunt resistance for each junction, and we found 
an astonishing result: In most of the ICFP, the effective 
shunt for each junction does not deviate more than 1% 
from Rq (Fig.25). 

The apparent success of a mean-field theory in this in- 
teracting system seems quite remarkable. Even more so 
is the fact that the RG-based calculation of the effective 
shunt on the critical manifold also yields with good accu- 
racy Rq. This is indirect evidence for an internal struc- 
ture in the coupled sine- Gordon theory that describes 
the two-junction system. If this approximation turns out 
to be generally valid, it may also give insight into the 
solution of coupled sine-Gordon models at intermediate 
couplings where no expansions can be done and little is 
analytically known. 

The relevance of our work goes well beyond the two 
junction system. A closely related class of systems con- 
sists of resistively shunted Josephson junction arrays of 
two or more dimensions (see e.g. Ref. 17). The RG equa- 
tions derived for the two-junction system^ were shortly 
thereafter also derived for a two-dimensional triangular 
array of resistively shunted Josephson junctions in the 
weak coupling regime: in Ref. 13, Tewari et al. show 
that the weakly coupled triangular array undergoes a 
superconductor-to-metal transition almost identical to 
the direct NOR-FSC transition of the two-jimction sys- 
tem (the RG equations for the two systems differ by one 
coefficient). As you may recall, this transition is in the 
interesting and novel ICFP region. Tewari et al. also 
considered the square lattice, in which the RG equation 
are third-order in the Josephson coupling, a difference 
which is unimportant near the ICFP. Hence, the work 
presented in this paper also verifies many of the results 



in Ref. 13. 

In the future, we plan to adapt the Monte Carlo 
algorithm such that it can be applied to sizable ar- 
rays of Josephson junctions. This would allow for 
the first time direct controlled investigation of the 
Kosterlitz-Thouless NOR-SC* (normal-superconducting) 
transition in one-dimensional resistively shunted and free 
Josephson-j unction arrays. Another interesting applica- 
tion for the algorithm would be Josephson-junction sys- 
tems coupled to quasi-particle dissipation^®. 
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APPENDIX A: CRITICAL EXPONENTS OF THE 
ICFP 

The behavior of the two-junction system in the ICFP 
region at intermediate energies is determined by the crit- 
ical properties of the unstable fixed point. Here we derive 
the critical exponents and principal directions of the RG 
in the weak and strong coupling regimes. 



1. Weak Coupling 



The weak coupling RG equations are given by 



dJ ^ R 



dl 



Rq 



(Al) 

(A2) 



where 



Near the unstable fixed point. 



r = -%v^, j; = ^u, (A4) 

VrR n 



we can linearize Eqs. (Al) and (A2) by writing 

J = J*+j, J+ = Jl+j+, (A5) 
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and thus obtain 



dl 



J 

j+ 




(A6) 



The cigcnvahics and eigenvectors of Eq. (A6) give the 
relevant direction and exponent, and also the irrelevant 
direction and its decay. For the relevant direction one 
finds 



and for the irrelevant direction 

A_ = - (^—w — \/w'^ + 8uvj^ 



1 



w 



2 \ H \ \j u \j u 
Near the unstable fixed point we expect 

aT-^+ +bT-^- , j+ ~ a+T-^+ +b+T-^ 



J 



(A7) 
,(A8) 

(A9) 
(AID) 

(All) 



where a and b are determined by the initial J*-"^ and . 
At very low T we are close to the J = J+ = fixed point 
and have 



(A12) 



If the system is superconducting, which means that ini- 
tially J(°) and are in the region in which they grow, 
then Eq. (All) still holds for intermediate temperatures, 

but as soon as the J's arc far from the ICFP, the system 
crosses over to the strong coupling regime, where 



as follows from Eq. (9). 



(A13) 



2. Strong coupling regime 

The RG equations for the strong coupling case are 

(A14) 



d( , RRq 



d^D . _ , tRq 2 



where 



RQ{R + r) 
2Rr + R^ ' 



w = 



i?2 + 2Rr 



- 1. 



(A15) 



(A16) 



As in the weak coupling case we linearize these equations 
near the unstable fixed point 



i?2 + 2rR 



RQVrR 



'uw, 



Cd 



i?2 + 2rR_ 

RqR 



(A17) 



writing 

C = C + ^, Cd = Cd + zd, 

and thus obtain 

d f z \ ( 1^' 



dl 



z 

ZD 




(A18) 



(A19) 



Eq. (A19) is essentially identical to Eq. (A6), for the 

weak coupling case, except for the change of u, w to w, u. 
Therefore all other results can be copied from the previ- 
ous section as well. The eigenvalues and eigenvectors of 
the matrix in Eq. (A19) are 



1 / r 

{z, Zd) = 1 1' 




and for the irrelevant direction 



{z, Zd) = 1, 




2 \ R 

Near the unstable fixed point we expect 

ZD ~ a+T-^+ + b+T- 



aT-^+ + bT- 



(A22) 
(A23) 

A_ 

(A24) 

just as for j and j+ in Eq. (All) for the weak coupling 
case. 



APPENDIX B: APPROXIMATE CALCULATION 
OF THE NOR-FSC PHASE BOUNDARY 

By using the result of App. A for the ICFP and the 
critical exponents and irrelevant directions, Eqs. (A22) 
and (A23), we can find an approximate implicit equation 
which can be solved for the critical i? as a function of 
r and Co (or Jo). The calculation is identical for the 
two limits; we will demonstrate it for the strong coupling 
limit. 

The calculation is based on the simple assumption that 
for r, R and tuned to the critical values, the critical RG 
fiow line between the point (^ = Co; ^£>=0) and the fixed 
point (^*, ^^) given in Eq. (27), is a straight line. This is 
illustrated in Fig. 29. The slope of the line representing 
the critical manifold is 



tan^ = \zd/z\ 




(Bl) 



where (z, zd) is defined in Eq. (A23), and with u and w 
(which are functions of r and R) given in Eq. (A16). The 
implicit equation for Rc is 



Cd* = (Co-r)tan0. 



(B2) 
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NOR 







FIG. 29: We calculate the critical i? as a function of r and 
by assuming that the flow from the initial state, (^o, 0) to the 
flxed point, (i^*, C,o) is a straight line (dashed-dotted line). 
This turns out to be a reliable approximation throughout the 
ICFP region. 



More explicitly it becomes 
i?2 + 2rR 



RqR 



Co 




APPENDIX C: SOLUTION OF THE CIRCUIT IN 
FIG. 21 

The effective lead-to-grain resistance is obtained by 
solving the circuit in Fig. 21 while replacing the Ohm- 
meter with a current source providing a current /. The 
following equations express zero potential drop along the 
two loops in the circuit: 

(/ + h)R+(I + h- h) r + haC + (/i - h) aCl = 0, 

(CI) 

hR + haC - [h - h)aCl = {I + h- h) ■ (C2) 
The solution for Ii and I2 is 

{2R + aC-) +R{R + aC+ a^) 



h = -I 
h ^ I- 



{R + aC) (2r + i? + aC^ + 2aCl) 
raC^ — RaCli 
{R + aC) {2r + R + + 2aC|,) ' 



(C3) 
(C4) 



and the potential drop on the junction connecting the 
grain B and the lead C is 



(C5) 



Substituting Eqs. (C3) and (C4) into Eq. (C5), and di- 
viding by /, we obtain the effective resistance measured 
in the Monte Carlo calculations: 



-Rcff — 



{R/a + C^) {2r + R + a^^ + 2aC,l) 



(C6) 
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